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Transverse Spin Diffusion 
in a Dilute Spin-Polarized Degenerate Fermi Gas 
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We re-examine the calculation of the transverse spin-diffusion coefficient in a 
dilute degenerate spin-polarized Fermi gas, for the case of s-wave scattering. The 
special feature of this limit is that the dependence of the spin diffusion coefficient on 
temperature and field can be calculated explicitly with no further approximations. 
This exact solution uncovers a novel intermediate behaviour between the high field 
spin-rotation dominated regime in which D± oc H~^, D^^ oc T~^, and the low-field 
isotropic, collision dominated regime with D± = oc T"^. In this intermediate 
regime, y oc T^^ but D± ^ Dy. We emphasize that the low-field crossover 
cannot be described within the relaxation time approximation. We also present an 
analytical calculation of the self-energy in the s-wave approximation for a dilute 
spin-polarized Fermi gas, at zero temperature. This emphasizes the failure of the 
conventional Fermi-liquid phase space arguments for processes involving spin flips. 
We close by reviewing the evidence for the existence of the intermediate regime in 
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experiments on weakly spin-polarized ^He and ^He — ^ He mixtures. 
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I. INTRODUCTION 



The unusual features of spin dynamics in spin-polarized quantum systems have been 
intensively studied since the pioneering paper of Leggett and Rice |l|] on spin diffusion 
in normal liquid ^He. The main effect arises from the observation that the presence of a 
molecular field (induced by the applied magnetic field) leads to an additional precession of the 
spin-current which in steady state acquires a component perpendicular to the magnetization 
gradient; through the continuity equation this results in an anomalous reactive component 
(damped spin-wave) to spin transport. This "spin-rotation" effect is also present in the case 
of spin-polarized Boltzmann gases |^. 

From a microscopic point of view, it was expected that a number of qualitatively novel 
phenomena might arise in the case of sufficiently high polarizations 0. A natural suggestion, 
made by Meyerovich was that low-temperature spin-diffusion becomes highly anisotropic 
for finite polarizations. More precisely, processes involving spin-flips make use of the phase- 
space volume between the two distinct Fermi surfaces (for "up" and "down" spins) leading 
to a finite scattering rate in the limit of T 0. This is in contrast with processes involving 
scattering in the vicinity of each of the Fermi surfaces which are subject to the phase-space 
restrictions of unpolarized Fermi liquids and are thus characterized by the conventional 
Fermi-liquid behaviour of scattering rates, oc Meyerovich's suggestion was recently 

supported by measurements of the transverse spin-diffusion coefficient in weakly polarized 
liquid '^He [^. Theoretically, D± was calculated in the dilute gas limit at T = 0J^; while 
for T 7^ the only available estimates are based on a variational solution of the Boltzmann 
equation [|^,|| , or on a phenomenological Mathiessen-type rule [|10 . 



In this paper we present the analytical calculation of the finite temperature behaviour of 
D± for a dilute Fermi gas in the s-wave approximation, by solving the appropriate kinetic 
equations exactly. While at H ^ T our results for D± agree with those of previous work 
I^J^, in the opposite case of T ^ if the exact solution of kinetic equation results in a 
qualitatively new behaviour of D±. We show, that the crossover between the high- field and 
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low-field behaviours is described by two dimensionless parameters rather than only one as 
implied in Refs. [[7|,p|,p^ : The first is the ratio H/T, which determines whether the phase 
space that is available for interparticle scattering is due to the field-induced splitting of the 
Fermi spheres {H ^ T) or to the temperature smearing of the Fermi distribution function 
(T ^ H). The diffusion coefficient is, roughly speaking, inversely proportional to the square 
of the available phase space volume; and hence in the two limits we obtain D± oc and 
D_L oc T~^, respectively. The second parameter, 

C = {H/T){apF)-\ej./T), (1) 

[a is the s-wave scattering length and ep and pp are Fermi energy and momentum in the 
absence of magnetic field) is the ratio of two frequency scales appearing in the kinetic 
equations for transverse spin diffusion, ( = u!i/u!2- The scale 0^2 is determined by the 
thermal term in the transverse relaxation rate, 002 ~ 1/tt ~ ma?T'^, where m is particle 
mass; while uji ~ appH is the molecular field contribution to the precession frequency of 
the transverse spin current. (Its origin is the spin-dependent exchange scattering amplitude 
of two particles interacting through a spin-independent interparticle potential Thus, 
the value of ( measures the number of oscillations of the spin current during a collision 
time. For C ^ 1 the "spin-rotation" effect strongly affects the spin diffusion process, which 
remains anisotropic even for H <^T: in this case, we obtain D^^ 0.89i^||. The crossover 
to isotropic behaviour occurs for ~ 1, so that D_l ~ D\\ only once ^ 1. 

We note that the crossover between isotropic and anisotropic behaviours of -D±,|j at 
~ 1 has already been discussed previously [|lT],|12|. However, the qualitative changes in 



the collision integral which occur at ~ T were overlooked. Similarly, the subtleties that 
arise in the solution of Boltzmann equation at C ~ 1 have not been noticed in Refs. [[7|,P,p!0 



which however provide the correct results for C ^ 1. In the present paper, we will give the 
complete description of field and temperature dependence of D_l for a dilute gas (within 
the s-wave approximation for the two-body scattering) which is valid for any value of H 
provided that the gas remains degenerate for both values of spin projection. 



where p| ^ are the Fermi momenta for spin-up and spin-down particles: 

N-^^l are the numbers of spin-up and spin-down particles per unit volume, My is the longi- 
tudinal magnetization per unit volume induced by the field H, and is the density of the 
gas. A short account of the present work was already published in Ref. [1^. 



We emphasize that throughout this paper we concentrate on the hydrodynamic regime 
defined by finite longitudinal and transverse relaxation times, r_Lj|, which, in principle, satisfy 
the condition -D±,||fc^T_|_ y <^ 1, where fc~l/L^O(Lis the size of the system) is the 
characteristic wavevector of magnetization fluctuations. The T ^ limit of our calculations 
leads to different physics from that of the strict T = 0, collisionless regime discussed by I. 
A. Fomin JI^. Indeed, one should interpret our T = results as describing fluctuations 
with sufficiently long wavelength, and at low but finite temperatures, T <^ H, such that 
the hydrodynamic condition is satisfied (note that and ry diverge at T 0). The 
interesting issue of the validity of the linearized Boltzmann equation approach in the context 
of spin-polarized Fermi gases at low temperatures, raised in Refs. |ll^JT5[| , is left for future 
investigations. Here, we limit ourselves to the linearized equations, which should be valid 
even for T H as long as the external probe used to drive the system out of equilibrium 
is weak enough and sufficiently slowly varying in space. 

The plan of the paper is as follows: In Section || we discuss the zero temperature limit 
of the kinetic equations derived more generally in Appendix B. To get a better feeling for 
the physical content of these equations we begin with an analysis of the relaxation time 
approximation. We then solve the T — equations exactly and recover the variational 
result of Ref. 0. This is followed in Section |T| with the discussion of the finite temperature 
equations in the limit of small polarization, H <^ ep, and the resulting exact solution for 
the diffusion coefficient. Since the results of Section |I| are valid everywhere in the region 



H ^ T our expressions are thus covering the entire range of field and temperature values in 
the degenerate regime, T <^ ep. This analysis uncovers the two crossovers already discussed 
above. In Sect. our results are compared with the available experimental data in '^He- 
^He mixtures [0,0,0,|T^ and in ^He 0. The derivation of the self-energies (which are used 
as an input into the Boltzmann equations) is relegated to Appendix A. The calculations are 
instructive as they uncover the simple phase-space origin of the finite value of D± at T ^ 0. 

Throughout we use the units with ks = fJ'O = ^ = ^, where ks is Boltzmann constant and 
/xo/2 is the magnetic moment of a particle; also we take the volume of the system, V = 1. 

II. TRANSVERSE SPIN DIFFUSION AT T ^ AND THE TRANSVERSE 

RELAXATION TIME 

We begin with the kinetic equations for the transverse (off-diagonal in spin space) com- 
ponent of the density matrix in the rotating frame, n|j(t,r, p) (see Appendix B): 

f d p d\ , ^ ^ 47ra . r , ^ d^p' , t u 47ra . 

xiN^-N^)nuit,f,p} = -z{[(s{?(eKp),p))* - {^f^ (e^ipj^p))*] nu(t,r,p-) + 

+ ^T^[t,r-^-^^,pynl-4)^ (4) 

where /x is the chemical potential and ^ denote the contributions in Eqns. 
( |A^| - [A14"D quadratic in the scattering length a; similarly, the transverse self energy, 



^U^^ 0'^'^ is given by Eqn. ( |B32| ). In the T — > limit, the Fermi distri- 

bution functions for spin-up and spin-down particles, nt and nt, can be approximated by 

It is easy to verify that the r. h. s. of Eqn. (^) does not mix different Legendre components 
of the quantity n||(t,r, p) as a function of momentum. Therefore, we will search for the 
solutions of the form 

nuit^r.p) = g{t,f,p) + /(t,f,p)cos?/', (5) 
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where is the angle formed by vector p and the x axis. It turns out that in the s-wave 
approximation the higher Legendre harmonics do not contribute to n-^i{t,f,p) (see below). 
Eqn. (D implies that, in the rotating frame, the transverse magnetization and spin current 
are given, respectively, by 

1 f d T) 1 

MZ{r,t) = Ml{f,t)-iMl{f,t) = -J n^^{t,f,p)—^ = i^io 9it,r,p)p^dp (6) 



and 



1 f d^p 
J±{r,t) = Jl{r,t) -iJl{f,t) = — / pn^{t,f,p) cos ip 



2m J --uv^.->^.--r ^27r)3 
f{t,f,p)p^dp, (7) 



127r^m Jo 

where the indices x and y refer to spin space. 

To find the transverse spin diffusion coefficient D± one has to calculate the value of the 
spin current J± induced by a slow spatial variation of the transverse magnetization. We will 
assume, without loss of generality, that VMj is uniform and is directed along the x axis: 

g{t,f,p) = {G, + G2-x)go{t,p), g^^I = ^ I 9o{t,p)p^dp, (8) 

where \ Gi + G2 ■ x 1 for any x inside the sample. 

The value of D± should then be found from the macroscopic constitutive relation ||l|,pO|, 
which can be written as 

where the quantity ^ is referred to as the spin-rotation parameter. The non-zero imaginary 
part of the denominator in (j^) reflects the fact that, due to the spin-rotation effect, the spin 
current is not parallel to the driving magnetization gradient. 

Substituting Eqns. and (|]) in the Boltzmann equation (H) and separating the Leg- 
endre harmonics, we obtain two coupled kinetic equations: 
(9 , , 1 p d ^ 2a 



d . , i p C „, ^ , -^a . /"^ , , 2 , / T 

T^9o{t,P> + ^ , ^ i^-^f{t,r,p) + % / go{t,q)q dq{nl-n, 

at Gi + G2 ■ X 3m ox vrm jo ^ 



-t{N^ - Ni)go{t,p) = -Irel[go] - ^Isr[go] (10) 



m 



and 



at m m 



-Irell[f] - ilsrl[f] 



where the hnear operators representing the "relaxation" and "spin-rotation" contributions 
to the colhsion integral are defined (for an arbitrary function W{p)) as 



Irei[Wm = [lmEf/(et(p-),p) - ImS;f (ei(p),p)] W{p) - 



m 



X 



(12) 



Isr[W{p)] = [ReSf)(eT(p),p) - ReJ:f^ ie^ip),p)] Wip) - 



m 



{p-p'f-Ak^ 



(2vr)e 



(13) 



where s = {p + p') /2, and 



(14) 



COS ll) COS ll) 

The simplest way of understanding the physical meaning of the kinetic equations ([IO|-[TI[) 
is to appeal to the relaxation time approximation for the collision integrals, Eq. ([l 



Irell[W{p)] = —W{p) , Isrl[W{p)] = -Q^^'^Wip) 



(15) 



where t±_ is the transverse relaxation time, while VP''^ can be identified as a second-order 
correction to the precession frequency of the spin-current (see below). Upon integration, 
Equations (p!0|-p!T|) then read |^ (see also Appendix B): 



r± 



(16) 
(17) 



where d = pjp^ and 



m 



'Sua 



m 



Mil 



(18) 



The main feature of Eqn. (^) is the precession of the spin current around the molecular field 
originating from the exchange contributions to the single-particle self energies. 
This does not appear in the equation for magnetization (Eqn. (0)) which only precesses at 
the bare Larmor frequency. In the simple case taken here, dM^/dt = dJ^/dx = and the 
solution to Eqn. ( [T7| ) reads: 



^lW = -Y^^Ml + Jl,,(0)exp 



T±. 



(19) 



Thus, the steady state current takes the form (^) with 



( Jj^^(O) represents the initial value of transient spin current which in our case is independent 
of X.) 

Having used the relaxation time approximation to illustrate the physical content of our 



kinetic equations (10^11) for the transverse density matrix we now proceed with the exact 
solution of these equations. 

For the purpose of computing the spin diffusion coefficient it is sufficient to consider the 
steady-state solution (^ with a time- and coordinate-independent functions f{t, r,p) = f{p) 
and go{t,p) = Qoip)- In this case, the first two terms on the 1. h. s. of Eqn. (|TUp vanish 
along with the first term on the 1. h. s. of ([TTl) 0. With the help of Eqns. (|A5 HA6D one can 
easily check that Eqn. (|10|) is then solved exactly by 

9oip)=nl-ni. (21) 

To proceed with the solution of the Boltzmann equation we note that, already on sym- 
metry grounds, we expect 

Ireiilf] ~ ma'H^f , hriif] ~ {N^'^afHf , (22) 



^The only coupling between the first and second Legendre components would involve df/dx which 
vanishes in our case in the steady state. This justifies the form assumed in Eqn. (||). 
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which also follows from explicit integration: after a rather cumbersome calculation we find 
that for any function W{p) that vanishes outside the region Pi < p < Pi (as required by 

/ Ireii[W{p)ydp=-^ F{x,d)W{xp^)dx, (23) 



where 



F{x, d) = ^x' - 2x^d' - Ix^d' + x^d' + ^d' + + 
35 5 35 5 

+—9{x^ + d^-l)- (8x' + d^- l){x^ + d^- 1)^/2 . (24) 
35 



Thus, the two terms on the r. h. s. of Eqn. (|TT| ) are small compared to the third term 



on the 1. h. s. which is of order N^^^aHf. Therefore, Eqn. ([TT|) may be solved iteratively 
by substituting 

^"'>--' 4.aW-iV.) '"^">'> 
on the r. h. s. . Upon integration, Eqn. (y) then leads to: 

D.^^^A^ I , (26) 

25m2 p3 _ p3 j^oo j^^^^ [pgQ{p)]p^dp ' 



15m Ireii\p9o{p)]p^dp 

Isrl[P9o{p)]P^dp 



(27) 



Irell\p9o{p)]p^dp ■ 

In turn, Eqns. (0), (pSf-P^), and (PB|-p7D can now be used to extract the final expressions 



for D± and ^, exact in leading order in N^/^a: 



where the diffusion time td is given by 

1 _ 16a2 r 21 



Td 1897rm 1 - d^ \ 64 ^ ' 

- -e{2d^ - 1) ■ (1 + ^d^){2d'' - ifl^] . (29) 
8 J 
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This proves that the result originally reported as variational in Ref. is in fact exact 
in the appropriate limit (we note that this result has been substantiated earlier through 
the numerical evaluation of the expressions in Ref. ||^). Equations (p8|-p9|) show that the 
transverse spin diffusion coefficient remains finite at T —>■ 0, as was already discussed above. 
In particular, in the limit of low polarizations, when 1 — d ^ mHpp^ -C 1, we find 

rn{H) = l^H-\ (30) 



i ma 



The detailed results for the spin diffusion time, td, are plotted in Fig. |l| (dashed line) 
together with variational bounds, r^, on the relaxation time, t±, which appears in the 



relaxation time approximation [21 



III. TRANSVERSE SPIN DIFFUSION AT FINITE TEMPERATURES AND 

SMALL POLARIZATIONS 



In this section, we extend our results to the case of finite temperatures in weakly polarized 
limit, within the degenerate regime, T ^ ep. The appropriate Boltzmann equation is derived 
in the Appendix B (see Eqn. ( [B33|) ). The only difference from the T ^ case of Eqns. 
- O) is in the form of collision terms, 



Irel[W{p)] = W{p) J A{p,p')^^ - J A{p',p)W{p 



(27r^3 



(2^ 



IsAWrn = W{p) j B{pJ'){nl, - 4,)^^ - (4-4) / B{pJ')W{p')^^ . (32) 



Here 



A{pJ') 



Tia^ ( (p2 + 



p 



m 



ni, + ni, - 



1 — exp 



2mT T 



and 



2aV 

B{p,p') ^ — 4j9f + Ip — p'\ In 

m \ 



\p — p' 


+ 2pp 


\p — p' 


-2pF 



(31) 



(33) 



(34) 
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(here pp = {2mepY/'^ is the Fermi momentum of the non-polarized system). Once the kinetic 
equations are solved, the values of and ^ should again be deduced from Eqns. 



We start with showing that the solution of Eqn. ([lOD with df /dx = again takes the 
form 



(35) 



The difference from zero-temperature case (see Eqn. (pi])) is that the Fermi functions on 
the r. h. s. of Eqn. (|35|) are now given by 

1 



n. 



hi 



exp 



+ 1 



(36) 



rather than by simple step functions. 

Since for the function gQ{p) given by Eqn. ( |55D the 1. h. s. of Eqn. (|iy) vanishes, and 



so does the operator /sr[5'o] on the r. h. s. (see Eqn. (|32D ) we have to prove only that the 
function (7o(p) solves the equation 



/rdb(p)]=0. 



Using Eqns. (|31|) and (^), we rewrite (^) in terms of dimensionless variables. 



(37) 



2mT T ' 2T 



(38) 



as 



g{(T){ri + a) x 



(39) 



where g{f]{p)) = g{p)- It is convenient to transform Eqn. 

g"{k) - {h^ + Ti'^)g{k) + 2Ti'^sec\i^Tik g{k) = 
sinh ^ 



da = 



- 1 . 

55|) into a differential equation 



27r II (7 



cosh Tik 



TT sin /c/i sinh irk — h cos fc/i cosh wk) 



(40) 



for the function 
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/CO fj 
e'^'^gir]) cosh j-diq, 

with the normahzation constant \\ g \\ given by 

/■OO 

II 9 11= / 9{'n)dV- 



Using a function goij)) of the form p5|), we obtain 



^0^^) = e^-l+i - e^+^^ + l ^^^^ 
and the corresponding Fourier transform, 

, ^ 271 sinh ^ cos kh , , 

Finally, verifying that goik) is a solution of Eqn. (^Oj) proves that the function (|35D satisfies 
Eqn. (p. 

We can then re-write Eqn. (|lTD in the steady state limit as 

G2VF ■ (4-4) - ^^^'Vb) = -Irell[f{p)] - llsrl[f{p)] • (43) 



Here vp is the Fermi velocity and O!-^^ = 2apFH/n (cf. Eqn. (|T8D). 

With the help of Eqns. (pi]) and (p3D we obtain the explicit expression for the integral 
of the relaxational term on the r. h. s. of Eqn. 



roo roc 

/ Ireii[W{p)]p'dp^p% F{T,H,p)W{p)dp, (44) 
Jo JO 



F{T,H,p) 



2m ^ 4 



(45) 



Here the function W{p) is assumed to be localized in the vicinity of the Fermi surfaces, 
where vp \ P — Pf |~ max(T, H). 

This yields the expected estimate Ireii[f{p)] ~ ma?f{p) max(if^,T^) for the relaxational 
term on the r. h. s. of Eqn. (^5]). In the case of spin-rotation term, the functions nt" 4 
Eqn. (^) carry a factor of if, and thus Isri\f{p)] ~ {ppdYH f{p) <^ Vl^^^ f{p). 
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The ratio of the precession term Q^^^ f{p) to the relaxational term is always large for 
H T, whereas for ^ T it is of the order of ( (see Eqn. (P) and can vary from C ^ 1 
(high-field region) to C ^ 1 (low-field region). On the other hand, the field and temperature 
dependence of relaxational term itself is determined by the ratio H/T. Below we will use two 
different methods to solve Eqn. (^3]) in two regions, C ^ 1 ^iid H/T <^ T/ep, assuming in 
both cases that H,T <^ ep. Owing to the presence of a large factor (app)'^ in the definition 
of (, the two regions overlap. Since the highly-polarized case H ^ T (with no requirement 
H <^ ep) is described by the T — > limit of Sect. we will thus cover the entire degenerate 
region {T <^ ep) for all values of H. 

High-Field Behaviour, 3> 1: In this case the relaxational term in (^) is small, /re/i[/] ^ 
i^^^-*/, and the r-approximation estimate for the spin- rotation parameter (see Eqn. (|20|)) 
yields ^ ^ 1, so that in the spin-space the current is almost perpendicular to the 
magnetization gradient. In this limit the solution of Equation (|43|), fnip), can be obtained 
iteratively in complete analogy with the zero-temperature case of Sect. |I|: 

fD{p) = -^^^2 (nt - n\;) - -pfypGa IreiiWp - n^p] - ^ (J(f))2 ^2 hriWp - 4] • (46) 
Upon integration this leads to 



J I = —I- 



where J±[W^(p)] = J W{p)p^dp/{12'K'^m) (see Eqn. (0)), and the functional l](^)[W^(p)] and 
T±[Vr(p)] are defined by 

poo roo 

/ Wip)p'dp / Isrl[Wip)]p'dp 

r±[W{p)] = — , n^'\W{p)] = . (48) 

/ Irell[W{p)]p''dp / W{p)p^dp 

Jo Jo 

t I 

The values of these functionals at the exact leading-order solution, nj^—n^^ represent respec- 
tively the second order correction to the spin current precession frequency and the "diffusion 
time" Tp,. Calculating the value of Tp, with the help of Eqns. (^^ -^5]) and performing explicit 
integration in we arrive at ||22||: 
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3 8ma^{H^ + An^T^) 
^ . + = jj^^j^f^ll + ^ap.d - 21.2)1. (50) 

The first thing to stress is that in the entire high-field region the value of D± in ( ^91) 



differs from that of the longitudinal spin diffusion coefficient p3 



2 

D|l = !^ ^ (4/87rma2T')C(-l/3) , (51) 
where C(— 1/3) ~ 0.843 is the Brooker-Sykes coefficient [24|. The two diffusion coefficients, 



Dj_ and D\\, are plotted in Fig. 0. Note that the spin diffusion remains anisotropic even at 
T ^ H (provided that C ^ 1), in which case we obtain from (^): 



in agreement with Ref. [|TT|. The difference between D± and at ^ 1 reflects the fact 
that in this regime the spin rotation effects dominate, as can be seen from the leading role of 
the precession term —i^l^^^f(p) in Eqn. (^31). Obviously, no such term enters the Boltzmann 
equation for the longitudinal spin diffusion, and the longitudinal spin current is by definition 
parallel to the gradient of My. Thus in spite of the identical temperature dependence of D± 
and the nature of the spin transport differs in the two cases. 

The crossover between and behaviours of D± occurs at T ~ i7. Eqn. (^Uj) 
shows that the diffusion rate r^^ obeys "Matthiessen's rule" , 

111 

— = — + -, (53) 

Td Th Tt 

where the quantities 

'" = ^^^^ = 32^^ " 

can be thought of as the relaxation times for field- and temperature-dominated collision 
processes. 
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We note that it has been commonly accepted 0,|,|T^] that the crossover of D_l at ~ T 



corresponds to transition from anisotropic to isotropic behaviour of spin diffusion. We will 
see that in reality this transition occurs at much lower fields, H ~ {apF)T'^ /ep- 
The dotted line in Fig. ^ corresponds to the quadratic fit 



nvl II 



which is obtained by substitution tt ry in Eqn. (|53D and is widely used to describe the 
temperature dependence of D±. One can see that while this fit provides correct order-of- 
magnitude estimate, it deviates from the exact result ( ^9] ) at T ~ iJ. 



Crossover and Low-Field Behaviour, H/T ^ T /tp: In this region, the spin-rotation term 
on the r. h. s. of the steady-state equation ( ^3] ) is small in comparison to both the relaxational 
and precession terms, 

Isrl[f{p)] « Irell[f{p)] , Isrl[f (p)] « ^^'^fip) , (56) 

and thus, it may be omitted in the leading-order calculation. Furthermore, one may also set 
H = in evaluating the collision term /ren- 
in terms of the reduced variables (pH]), Eqn. (^31) then becomes 



vp ^ f I I \ 2o?m ^ / 2 , 2 

4a^m r"^ , . „, . / 1 1 



, (7] + a)f(a)[ + — ]da, (57) 

where f{r){p)) = f{p) and h <^ 1. Since the 1. h. s. of this equation is even in r] and the 
integral term on the r. h. s. has the parity of the function /(r/), the inhomogeneous equation 
(^) should have an even solution /(r^) = /(— 77). In finding this solution, we will make use 
of the methods of Ref. and begin by transforming (pTD into the differential equation 



^»(,) _ ..,^^(,) - ?.ch^.. Fi,k). ^ (58) 

for the function 

F{k) = r e*''/(^) cosh ^ dr] , F{k) = F{-k). 
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Here 

The solution of Eqn. (|58|) can be written as an expansion 

oo 

F{k) = Y.Fr,<t>n{k) (59) 

n=0 

in the basis provided by the functions 

Uk) = {l- X^V'^CT^"" (X) , X = tanh vrA; , (60) 

where Cn^^(x) are Gegenbauer polynomials with a complex index 7 + 1/2. The functions 
0„(A;) satisfy the homogeneous equations 

^"{k) - 7r V0(^) - (7 + Ti) (7 + n + l)sech27rfc <\){k) = (61) 

and are orthogonal to each other, 

[ (j)nik)(f)ra{k)sedi^TTkdk = ^(^7 + ^ + ^) (62) 

n! (n + 7 + I) [r(7 + I) 

Since (pn^—k) = (— 1)"0„(A;), only the functions 0n(fc) with even n contribute to the expan- 
sion (^) of an even function F{k). Substituting the expansion 

— — - ~ — - = secli^Tck gn<Pn{k) , (63) 

3 coshTTK 3 coshTr/c 

where the summation again has to be performed over even n only, and 

7r5/2(n + 7 + i)r(7 + |) 



^" 3cosf r(i - f)r(i + 7 + f)r(i + f + f)r(i - i - 

into the r. h. s. of Eqn. (|58D we finally obtain 

3vph 



(64) 



^G2 , , , J^:, , (65) 



27ia^mT^ ^ | + (7 + n)(7 + n + 1) 

for even n. As was already mentioned above, F„ = = for odd n. 
The transverse spin current (0) can now be evaluated as 
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00 



an ^ I ^4^dk = (67) 
_ 7rr(7 + l)r(27 + n + l) 

" n\ cos ^ r(27 + i)r(i - f )r(7 + f + i)r(| + § + i)r(i - f - i) 

for even n, and = for odd n. Finally, Eqns. (|^ - and (^) yield the expression 



00 



3HG 



a. 



2 „=o 



~ 4^;;^^^J + (7 + n)(7 + n + l) ' ^ 
from which one can extract D± and spin-rotation parameter ^ (see Fig. |^). In the high-field 
limit C ^ 1) Eqn. (|68D gives 

~ 0.890D|| , e = -^C > 1 (69) 



in agreement with Eqns. ( p2D and (|50D, while in the low-field case of C ^ 1 we obtain the 
isotropic result, 

D± = Dii, ^ ^ 0.139C < 1 . (70) 

Note that while the value of D± in (|70D can be obtained from the r-approximation formula 
(^) with r_|_ = Ty (see Eqn. (|5TD ), spin-rotation parameter ^ is not given by 

The crossover between Eqns. (|69|) and (^) at C ~ 1 cannot be properly accounted for 
within the relaxation-time approximation (Eqns. (pO|)) because the relaxation operator Ireii 
on the r. h. s. of Eqn. (|43| ) (see also Eqns. (^^-|45|)) does not depend on the magnetic 
field for H <^T . This is in contrast with the interpretation of the low field crossover as a 
result of the competition between two different terms in the relaxational part of the collision 
operator As illustrated in Fig. ^, what changes at ~ 1 is the form of the solution, 
f{p), to the steady-state equation (|43|). It is this change in f{p) which is not accounted for 
within the relaxation time approximation, hence the absence of the second crossover in that 
approach [|T0[] . 



Although quantitatively our considerations apply only to the degenerate case (T -C ej?), 
there is no reason to doubt that the crossover at ^ ~ 1 should also be present in the 
case of ~ {apF)eF and T ~ ep, in the presence of higher partial waves, or in a system 
with stronger interactions. Since the separation of scales of the two parameters is provided 
by the factor [ep /T) / [app) ^ 1 in Equation (|I|), in both of these cases the intermediate 
region between the two crossovers should shrink. In the next section, we will discuss the 
experimental results that support this point of view. 



IV. DISCUSSION OF EXPERIMENTAL RESULTS 

Some comments are in order concerning the possible relevance of our findings to experi- 
ment. Even though, strictly speaking, our analysis does not apply to the strongly interacting 
case, it is worth noting that the available experimental data in weakly polarized ^He 
deviate systematically from the simple theoretical fit which uses a single adjustable param- 
eter {Ta in pi) to cover the entire temperature range including both Q > 1 and < 1. 
Much better agreement is obtained by restricting the fit to the ~ 1 region with an overall 
prefactor smaller than the one implied by fitting to the value of D\\ in the low field, high 
temperature regime (in re-analyzing the data we use the value A 7.7 x lO^sec/cm^K^ 
instead of 5.8 x lO^sec/cm^K^ for the inverse overall prefactor in Eqn. (2) of Ref. 0; we also 
take the value 12.5 mK instead of 16.4 mK for the "anisotropy temperature" Tg). This is 
consistent with our picture, with Dx_ < D\\ for ^ ~ 1 (i. e. ^ ~ 1) and T ^ H. In addition, 
although it appears that the region between the two crossovers cannot be clearly identified 
- most likely due to large Fermi liquid renormalization effects - the isotropic limit is indeed 
reached in the regime ^ < 1 0. We also note the systematic deviations from the quadratic 
fit of Ref. ||10[ visible in the measurements of D± in a 6.4% '^He -^He mixture |T6[ in the 



intermediate region between T Ta and | ^ |~ 1. However, the interactions in this system 
are still too strong for our theory to be applicable. The fact that in both cases of Refs. |0 
and [TO] the intermediate region appears to be rather narrow is in line with our expectations 
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(see Sect. In principle, our calculations should be more relevant to the measurements 
in dilute ^He -^He mixtures. Although in the available data (for .18% '^He) the crossover 
to the isotropic limit occurs for ^ ~ 1 with H <^ T, the temperature is not sufficiently far 
below ep and, moreover, the polarization is somewhat high, ~ 25%. Nevertheless, for > 1 
D± (X D\\ with the ratio D^/ D\\ slightly less than unity |0,|T3. Also, the measured T depen- 
dence of the "spin-rotation" parameter, ^, near the crossover to the isotropic ("low-field") 
limit is quahtatively consistent with our results in both the data of Reference [jl^ and those 
obtained in the degenerate regime of more concentrated solutions (2.6% ^He) with lower 
polarization (~ 2%). In both situations, the crossover to the isotropic regime can be clearly 
distinguished. However, in the former case the polarization was again rather high (~ 25%) 
and the "low- field" crossover itself occurs beyond the degenerate limit. Also, there is a large 
discrepancy in the magnitude of the shift of ^T^ in the latter case (see Fig. ^ which can 
be attributed to Fermi liquid renormalizations anticipated in high concentration solutions. 
The deviations of ^ at ^ ~ 1 from the predictions of Ref. are also seen in the data for 
the dilute (.05% and .1% ^He ) mixtures of Ref. |TB|, in which case due to the high polar- 
ization ^ ~ 1 corresponds to T ~ e^?. Again, our results do not apply quantitatively to this 
situation. In addition, the data points are somewhat too scattered for a clear identification 
of the crossover. To sharpen the identification of two crossovers the data of reference 
should be extended to lower temperatures (to study the H/T ~ 1 behaviour). Quantitative 
comparisons could be made only in the more dilute case of references ||I7|JT6|1 where lower 
field and lower temperature experiments should be performed. 
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APPENDIX A: ON THE DYNAMICS OF QUASIPARTICLES IN A 
SPIN-POLARIZED FERMI GAS 

In this Appendix, we derive the imaginary parts of quasiparticle self-energies at T = 
which are used as input into the kinetic equations. 

We start with a spin-polarized gas of neutral fermions interacting via a repulsive potential 
U{\ f ~ f' I). In the dilute limit, N^^'^r^ <^ 1, where tq is the range of the potential, 
the leading contribution comes from the s-wave scattering between spin-up and spin-down 
particles. Then the natural small parameter is provided by N^^^a, where a is the s-wave 
scattering length. We will keep terms up to the second order in N^^'^a; within this accuracy, 
the scattering is fully described by a two-body T-matrix which obeys the Bethe-Salpeter 
equation: 

T{pi,p2]{so,s)) = -U{pi-p2) + / -nrx?- — — — — — TTnTfT^ ' 

J [l-Ky 2sq - s ^/m + 2^1- k^/m + iM^^ (k) ■ 

where 

f^^^(k) = 1 - nl .-ni . , (A2) 

^ ^ s-\-k s—k ' ^ ^ 

T I 

and ni are the Fermi distribution functions for spin-up and spin-down particles respec- 
tively, and (so, s) is the 4- momentum of the center of mass. The momenta of the incoming 
(outgoing) particles are s±pi, (s±p2)- Equation ( [Al[ ) is solved by 

^ ^ Ana (ATiay f Af^^{k) 

T{pi,p2]{so,s)) 



-V 



m \ m 
m 



pi - k"^ 



2so - s Vm + 2/i - k^/m + iAf^^{k) ■ 
" (A3) 



(27r)3 ' 

which is a straightforward generalization of a well-known result [PB| for the case of a spin 
polarized system. 
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To obtain the self energies of spin-up (T,^^{u,p)) and spin-down {T.^^Ujp)) particles, 
one has to perform the integration of the T-matrix multiplied by the single-particle Green's 
function, as shown in Fig. ^ We will first consider the on-shell self-energies Ti^^{uj = e^{'p),'p) 
and = e|(p),p), where 



are the single-particle energies. To second order in aN^^^ we obtain 



Ana 



m 



Ni + 



(Ana) 
2m 



lip- P'f + imM^^{k) ■ 



24, 



47ra 



(47ra) 
2m 



U^i(k){U^^{k) + 2nl,-l 



24, 



Up-p')2-k2 + imArn{k)-0 \{p-p' 



A;2 



d?k d?p' 
(27r)6 



d^k d?p' 
(27r)6 



Here the center-of-mass momentum s that enters the definition of N'^^{k){see Eqn. 
equal to {p + 'p')/2. 

Separating real and imaginary parts in Eqn. (K^) yields 



Ajrn 

ReSn(p) = ^iV^ + 



ina] 



m 



m 



V 



s+k s—k P s+k y s—k 

{p — p'Y — 4fc2 



d^k d^p' 
{2Txf 



and 



ImS||(p) = vr 



(47ra) 



m 



i 

s—k 



X 



2 {p~P'y\ d^kd^p' 



(A4) 



(A5) 



(A6) 



IS 



(A7) 



(A8) 



4 ; (27r)6 ' 

similar equations with t*^i hold for ReS|| and ImS||. In the non-polarized case, when 



n^ = , the expressions (|A5| 



reduce to the familiar formulae for the self energy p6| , p7 



At T = 0, it is possible to perform the integration in Eqn. ( [A8| ) (and in the similar 
expression for ImS^) explicitly. After a cumbersome calculation we obtain the following 
results for the imaginary parts of self energies: 
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At g > 



ImS||(g) 



r 

~ 15^ pg^^ - 3pt' + Pi' - WPi'+ 
+ 2e{p,'+p,'-q').{p,'+p,'-q'f'' 



(A9) 



(AlO) 



In the intermediate region Pi < q < p^ 



ISvrgm 



2\5/2' 



ISvrgm 

And, finally, in the innermost region q < Pi 



2g5 - 5q\^ + ?>p{> 



(All) 
(A12) 



ImS||(g) 



GOvrgm 



2\2 



Equations ( [A9p and ( |A14| ) coincide with those reported in Ref. 
[A13|) are presented here for the first time. 

The singularity at = Pt^ + Pi^? which is present at Eqns. 



(A13) 

(AM) 

The expressions ( [A1(J| - 



9 HA10|) corresponds to 

a process in which an incoming particle excites a particle of opposite spin from within the 
Fermi sphere. At g > (p|^ +Pi^)^^^5 the incoming particle has enough energy to excite 
any particle of opposite spin. Geometrically, this singularity corresponds to the exterior 
tangency of the shifted Fermi spheres, and it is already present in the unpolarized case . 

In the spin-polarized case there is one more singularity at g^ = — Pj,^ (see Eqns. 
( |A11| ) and ( |A13| )), associated with the interior tangency between the two Fermi spheres. In 
this case, the corresponding virtual process is a recombination of the incoming spin-up hole, 
followed by an excitation of spin-down particle; at g < (p|^ — pj,^)^^^, the hole has enough 
energy to excite any of the spin-down particles from within the spin-down Fermi sphere. 



23 



The functions ImS||(g) and ImS||(g) change their signs at q = and q = Pi respectively. 
The new feature of the above results, relevant for the transverse spin dynamics, is that, due 
to the availability of the scattering phase space between the two Fermi surfaces, the lifetimes 
of spin-down particles and spin-up holes in the intermediate region Pi < q < Pi are finite. 

The transverse spin processes are described by the Wigner transform of the transverse 
component of the density matrix, n||(t, r, g), which obeys the Boltzmann equation of a 
peculiar type (see Sect. p|-pD. In the appropriate linearized approximation sufficient 
for the calculation of transport coefficients the transverse components decouple from the 
longitudinal ones, enabling one to discuss the transverse and longitudinal spin dynamics 
separately ||2^. By rotational symmetry, the transverse deviations from the equilibrium 



(uniform) polarization must be small at any point in p space, i. e. 

nnit,r,q} <^^ni-ni. (A15) 

This also respects the fact that the eigenvalues of the density matrix are between zero and 
unity. 

Thus, in a non-equilibrium state the value of n-^i [t, r , q) has to differ from zero over a 
finite interval of values of q within the intermediate region p^ < q < pp The fact that the 
quasiparticle relaxation rates ( |A11| - [AT^ ) also differ from zero within this region results in 



the finite relaxation time t± of the transverse spin current even at T ^ (see Sect. IH). 

We note that from the viewpoint of longitudinal spin dynamics the system behaves 
as a conventional Fermi liquid: the non-equilibrium terms in the diagonal components of 
density matrix are confined to the respective Fermi spheres, where the imaginary parts of 
the quasiparticle spectrum vanish. As expected p8[, the imaginary part of the off-shell self 
energy at the Fermi surface behaves as 

2 

ImS||(e|(p|) +iu,pi) ^ — signu; ■ u;^ . (A16) 
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APPENDIX B: DERIVATION OF BOLTZMANN EQUATIONS 

In this Appendix, we derive the kinetic equation for transverse spin processes using the 
Keldysh approach We will consider time dependent Green's functions of the form: 

^Gta2p = {'^iMp) , = -(^2/3^1a) , (Bl) 

which satisfy the identity 

iG-,-,p + - ^Gt-^p - ^G-,^2p = . (B2) 

Here, the indices 1,2 stand for the two space-time arguments Xj = {ti,fi), T and T are 
the operators of chronological and anti-chronological ordering respectively, and the average 
(■ ■ ■) is taken over the given non-equilibrium state. 

In the absence of the interaction, the functions (|B1|) obey the following equations of 
motion: 

GoiajGf^2t = Tf6al3S{Xi - X2) , (B3) 

* Gf^^i = r:%p6{X, - X2) , (B4) 

where the indices a, b accept values + or — and r"'' is a Pauli matrix: = — t++ = 1, 
T^^ = = 0. Here and below we presume that the sum should be taken over the repeated 
Latin or Greek (i. e. spin) indices. The operator Goiajs is defined as 

The Green's functions of the interacting system satisfy the Dyson equation 

or, equivalent ly, 

Gab ^{0)afe / ^ac s;^cd /^{0)<ib j4 -y- j4 -y tT}'7\ 

la2/3 — *-^la2/3 + J <-^la37^374A*-^4A2/3 " A3 A4 [tii ) 
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where S"^2/3 the proper self-energies. 

Our aim is to examine the equation of motion for G^^^^i which is the off-diagonal com- 
ponent of the density matrix in spin-space which is ultimately related to the transverse 
magnetization and associated spin current (see Eqns. (H - 0))- To derive this equation 
we act on Eqns. ( |B6D and ( P7| ) with the operators Gq^i , and (60211) , respectively, and 



subtract the resulting equations from one another to obtain, 

'-^im = ~ J ^lT3Q^3a2i + ^lT3a^3Q2i d.'^^S , (B8) 



d 1 

-i— V,-.V,- - H 

ot m 



where X and X' are defined by 

X = ^(Xi + X2), X' = Xi-X2. (B9) 

We will be interested in departures from equilibrium which are (i) sufficiently small so 
that the equations can be linearized and (ii) slowly varying as functions of X so that only 
leading terms in the gradient expansion are required. The first approximation implies that 
only one of the internal lines in any diagram contributing to S^2j^ or contain the 

transverse Green's function G'^^^^; by symmetry, the departures of the diagonal functions 
Gi^2'\y ^ii2i) ^iT2T ^ii2i from their equilibrium form are already quadratic in G3141 
resulting in the decoupling of the transverse and longitudinal Boltzmann equations in the 
linearized limit. The second approximation requires that, in the rotating frame, the char- 
acteristic scales of variation of G'^'^^ as a function of X should be large compared to all 
microscopic scales in the system: 

1 TTl 1 1 

AL > — , At > — , AL > , At > — . (BIO) 

Pi Pi Pi- Pi H 

The separation of scales discussed above is implemented as usual by first transforming 
all quantities into the rotating frame and then performing the gradient expansion in the 
Wigner transform. The latter is defined as the Fourier transform with respect to the relative 
coordinate, X' (here F{Xi,X2) is an arbitrary function): 
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(2^ 



F{X, P) = J e'^-^'F(X + X - ^X') d^X' . 
P = {uj,p), P-X' = ujt' ~pr' . 



(Bll) 



In transforming into the rotating frame all transverse self energies and Green's functions 
acquire Larmor precession phase factors, 



/^ab iHtf~iab spab iHtspab 



(B12) 



For example, for the mixed Fourier transform (only time arguments are shown) of a term 



/C containing 



^(^1,^2) — j ^{tl — — ti iti — t2 )G'3|4| dt^dti 



one obtains. 



1 



J^(uj--H,uj~uj\uj + -H) G?t (cj', t) 



duo' 



(B13) 



For the purposes of deriving the kinetic equations, the rotating frame function , ^(cc;, t), will 
then be treated as a slowly varying in the center-of-mass time variable, t. 

With the help of Eqns. (pil| - pT3D the kinetic equation ( p8|) can be rewritten as: 

( d p d 



dt '''' mdfj ^ ^ ' ^ 



- Sf^"(P_)G!f+(X,P) + 
+Sff (X, P)qi+(P+) + Gff (P_)S'^+(X, P) + Gff (X, P)SJ+(Ph- 



(B14) 



where P± = {u ± \H , p). Note that in Eqn. ( |B14| ), the dependence on X only appears in 
the transverse Green functions and self-energies as expected in the linearized limit. 

The various self-energies on the right hand side of Eqn. ( |B14| ) are related to the T- 
matrices as shown in Fig. since both ends of the interaction line bear the same sign {+ 
or — ), there are only four kinds of the T-matrices involved (T , T^+, and see 
Fig. 0). All T-matrices are determined by Bethe-Salpeter equations which, to lowest order 
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in the transverse Green functions, conserve the spin of each of the incoming particles. It is 
convenient to rewrite the corresponding integral equations by separating out the T-matrices 
involving independent propagation on each of the Keldysh contours, Tq and Tg*'^; the full 
T"*-matrix is then obtained by connecting factors of Tq and Tq'^ by the appropriate G^^ 
or insertions (see Fig. ^). 

Thus Tq and Tg^^ are determined by the Bethe-Salpeter equations with no G or 
G^^ intermediate lines. It is easy to see that, for T = 0, 

To~{pi,P2; (so, s)) = T(pi,p2; (so, s}) , (B15) 

where T(pi,p2; (so, s)) is the usual equilibrium T-matrix calculated in Appendix A (see Eqn. 
I)); while at T 7^ we have: 

4™ Mvra^ ^ 



^ Aira MvraX 

^0 {PuP2;{so,s}} = X 

m \ m J 



X V 



2so — s^/m + 2yU — /m p^ — /c^ 



+ (B16) 



(27r 



i3 



(see Equation (|A^)). Similarly, 



(2vr) 



To'^iPi,P2] (sq, s)) = - (Tq (pi,P2; (so, s))) . (B17) 

We note that the ( — h) insertions cannot modify the physical scattering length! and 
thus, as can be seen from Fig. ^ b, to second order in (N^^^a), the T-matrices T and T~^~^ 
are determined by the first term, 

T~~{puP2; ; (so,s)) = Tq-{pi,P2; {so,s)) , 
T++(pi,p2; (so,s)) = - (Tq'{pi,P2; (so,s))) . 



^Formally this follows from the fact that the integrals over energy and momentum in a "bubble' 
composed of two (— +) Green's functions are convergent. 
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Similarly, to second order in the scattering length, T~~^ is determined by the first diagram 
in Fig. H b, proportional to an integral of Tq Gj^Gj-^^Tf^^ . 

Finally, the self-energies required in the derivation of the kinetic equations can now be 
calculated from Fig. |^. For example, the real and imaginary parts of the on-shell self- 
energies (Fig. 1^) S^|^(e|(p), j)) are given by Eqns. ( |A7D and ( |A8D (see Eqn. ( |Bl(j| )); similar 
equations hold for S^^(e|(p),p). Moreover, 

S++(P) = - (Sfr(P))* , ^ttiP) = - (S,y (P))* . (B18) 



We are now in the position to reexpress ( |B14| ) as a kinetic equation for the transverse 



component of the density matrix (in the rotating frame). 



nnit,r,P) = -ij G7+(X,P)— . (B19) 

We first note that the expansion of the off-diagonal self energies S^^^ in powers of N^^^a 
starts with the terms proportional to (N^^^a)'^ (see Figs. In the diagonal self-energies, 

and Sq^, we separate the first-order contributions, 

(P) = .^N^ , (B20) 
m 





{P) = 


m 


(X, P) = 




m J 



and move the corresponding terms in Eqn. ( [B14| ). to the 1. h. s.. Expanding the Green's 
functions Gj^^ and G^i^ with the help of Dyson's equation ( [B6D and Eqn. ( [B20[ ), we find, 
to linear order in n-^i{t,f,p): 

S;y-(X, P)G^7(P+) + Gf +(P_)s(y++(X, P) = (B21) 



Here we also used the following formulae for the non-interacting longitudinal Green's func- 
tions: 



Gi?-(P_) = - (g{?++(P_))* = V ^— + Z7r(24- 1)5 
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2m 



(B22) 



G 



(0)- 

TT 



'{P-) = 27iin^J{uj 



2m 



and similarly for the J,| components, in which case P_ is replaced by P+. Keeping the terms 
of up to the second power of appi we thus re- write Eqn. ( |B14| ) as 



i6 {u - e{p) + fi) {nl - ni) I n^i {t, f, p') -j^^ + 



(2vr) 



+S{;)--(P_)G7+(X, P) + +(^, P)^T^iP, 



sS^""fP_)G'J°^"+(X,P)+ (B23) 



-Sjf-"(X, P)G'[°)"+(P+) + G'}';)-''(P„)S{|'"+(X, P) + ^{"^""(X, P)s[f +(P+) 



where to the leading order in interaction we can replace the transverse Green's functions on 
the r. h. s. with 



Gf^~''{X,P) = Gf/+"(X,P) = g{°/++(X,P) = Gf^''{X,P) 
= '2nin^i{t,r,p)5{uj - ^ + /^) ■ 



(B24) 



Finally, using also Eqns. ( P20|) and (P22| ), we find upon integrating Eqn. (|B23| ) over a;: 



47ra 



m 



■i(X| - Ni)n^l{t,r,p) = -Id[n^i{t,r,p}] - Ie[n^i{t,f,p)] . 



(B25) 



Above, Id[n^(t,r,p)] and /e[n||(t, r, p)] are the direct and exchange contributions to the 
collision integral, respectively: 



Id[n^l{t,r,p)] =i{s^+(e|(p),:p) + sj|^ {e^{p),p) + 



(B26) 



and 



h[n^^{t,r,p)] = 2 ^ - nt- 1) + 



2m 

u 



2m 



2m 



it 



(B27) 
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where the superscript (2) again imphes that only the quadratic terms in a should be included 
in the corresponding self-energies (we have used the notation ( |A4| ) ) . 

Equations ( |B25| - |B2^ ) are valid at any temperature; at T — ^ 0, the kinetic equation can 
be simplified further by assuming that 

nnit^r,p) = 0, p<pi orp>p|, (B28) 

This ansatz is consistent with the linearized equation: to see this, first note that Eqns. 
( [B22D , (|B24| ), and (|B28| ) together with energy conservation in the corresponding diagrams 



of Figs. 0, H & imply 



2 



'2m 

From the transverse component of the self-energy identity p9[, 



^la2/3 + ^la2/3 + ^la23 + ^la2d " (B30) 



it then follows that the exchange contribution to the collision integral ( |B27|) indeed vanishes 



outside of the intermediate region, Pi < p < (the direct term already vanishes by (P28|) ). 
Similarly, in the case of longitudinal self energies, energy conservation yields: 

SfT~^(eT(p),P) = 0, P>Pr, S^+(e|(p),p) = 0, p>Pi; 

Stt~(^t(p)>p) = 0, p<pt; ^tr{^i{p),p) = p<Pi 



which, with the help of ( p30D and ( [BTSD lead to: 



Stt^(^t(p)'P) = -2«ImE^^(e|(p),p) , p < p^ ; 

^lt(^l(P)^P) = -22lmE^(e|(p),p) , p<Pi. (B31) 

Equations (|B29HB31|) then enable us to re-write the Boltzmann equation at T — >^ in the 
form M). 
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The finite temperature Boltzmann equation is derived by substituting the exphcit ex- 
pressions, 



-TCI 



r 



s+k s—k 



(B32) 



(2vr) 



m 



s— A; s+fc 



{4,-i)s[\ip-pr-k') 



2^ d^k d^p' 
d^k d^p' 



(2vr) 



into ( P25| - p27|) (here and below s = {p + p')/2) to obtain the equation vahd at any temper- 
ature and magnetic field in the dilute s-wave limit: 



X 



Ana 
m 

Ana 
m 



- B{p,p') 
-B{p',p) 



1 d^p' 
{nl - nl) - A{p', p) \ (t, f, p') - 



(2vr) 



(B33) 



Here 

A{p,p' 
and 



B{p,p') 



m 



(B35) 



{p-p'f-Ak'^ (27r)3 

We note that Eqns. ( |B33| - [B35| ) are in complete agreement with the appropriate limit (see 
inequalities ( |B10| ) above) of the kinetic equation of Ref. [0^. Upon integration of Eqn. 



^The origins of the difference in numerical factors between Eqns. ( [B33| - B35| ) and the results of 
Ref. pH] (which contains another derivation of Boltzmann equation within Keldysh formalism) are 



unknown to us. 
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( P33|) over the momentum space the two terms on the r. h. s. ("direct" and "exchange" 
terms) cancel each other, and, according to Eqns. we are left as expected with the 

continuity equation (0). 

Elementary integration in Eqn. (|B34|) leads to the following expression for A{j),p'): 



A{p,p') _ 2 1 
2TTa'^T s 1 — exp(2?i;) 
1 + exp{su + w + h) 



In 



1 + exp(sM — w + h) 
1 + exp(— SM — w + h) 



-In 



1 + exp{—su + w + h) 



+ 



(B36) 



1 + exp(— — w + h) 
where we use the following notation: 



J 1 + exp{su — w + h) 1 + exp{su — w — h) 



1 + exp{—su — w — h) 



u 



p — p 
2mT 



w 



ijp' + p'^) — 4m/i 
4mT 



h 



H 
2T 



In the limit of low temperatures and small polarizations, when T,H ^ /i, only the 
electrons from the vicinity of the Fermi spheres may participate in the scattering. Therefore, 
in the expression ( [B36| ) we may assume that p ^ p' ^ (2m/i)^/^ and w ~ max(T, iJ)/r. 
Moreover, it is easy to check that the contribution of backward and forward scattering 
regions (where the product su may be arbitrarily small) is negligible, so one can also assume 
that suT ^ max(T, H). 

With this assumptions, we obtain Eqn. (|33D for A{p,p'). The quantity B{p,p') at 



T, <^ /i is given by Eqn. (0). 
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FIGURES 

FIG. 1. The upper and lower bounds of the relaxation time functional t_l as functions of 
d = Pi/pi- Solid lines correspond to the functions a^p-^^T^{d)/m . The dashed line represents the 
diffusion time (see Eqn. (|29|)). 



FIG. 2. Spin diffusion coefficient in the high-field region H ^ T. The transverse spin dif- 



fusion coefficient (Eqn. (^9|)), the longitudinal coefficient D\\, and the simple fit ||;1C| 
(l/nf^ + l/DiiY^ (Eqn. (H)) for D_l (of^ = 3TTv'j,/8ma'^H^ is the limiting value of D_l at 
T —f 0) are represented by the solid, dashed and dotted lines, respectively. 



FIG. 3. The crossover at T ~ Tc = Hep / app (i. e. C ~ !)• Fig. y a shows the ratio D_\_/D\\. 
The spin-rotation parameter ^ is plotted in Fig. ^ h. 



FIG. 4. The functions f{p) that solve the steady-state equation (14 3| ) in the isotropic limit 
Q <^ 1 (solid line) and in the intermediate regime, C 1 ^'^d H <^ T (dashed line) at the same 
temperature T. 

FIG. 5. The low- field behaviour of the spin-rotation parameter ^ in the concentrated ^He-^He 
mixture. The points represent the experimental data of Reference the dashed line - the 
fit for the behaviour of ^ below the crossover, and the solid line is our theoretical result for the 
s-wave approximation (with the values of T and ^ scaled by appropriate factors). Note, that the 
latter is not a straight line. 

FIG. 6. The diagrams for the self-energies E||(p) and T,i^^{p). 
FIG. 7. The T-matrix T"'^{pi,p2; {sq, s)) and the diagrams for the self energies S^i, S?!* and 



S|J in the Keldysh technique; the indices a and b may take the values -|- or — . 
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FIG. 8. (a). Bethe - Salpeter equations for the auxiliary T-matrices Tq and Tq~~^. (b). The 
diagrammatic series for the T-matrices T and r~+. The series for r++ and T"* can be drawn 
in a similar way. Within the second order in N^^^a only the first terms in these diagrammatic 
series should actually be taken into account. 
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